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ABSTRACT 


Reentry  vehicle  roll  resonance  due  to  small  mass  and  aerodynamic 
asymmetries  is  described  analytically  by-  a  quasi-steady  solution  of  the 
three -degree -of-freedom  moment  equations  of  motion.  The  quasi-steady 
analysis  neglects  oscillations  in  the  vehicle  orientation  and  angular  rates 
about  average  (quasi-steady)  values  of  these  parameters,  which  change 
slowly  with  time  relative  to  the  oscillations.  The  study  treats  different 
regimes  of  reentry  vehicle  roll  resonance  that  have  been  identified  through 
computer  solutions  of  the  general  equations  of  motion.  These  include  high 
altitude  roll  lockin  and  breakout  due  to  a  single  mass  asymmetry,  inter¬ 
mediate  and  low  altitude  sub  resonance  with  a  single  asymmetry,  and  spinup 
and  lockin  to  low  altitude  resonance  from  a  compound  asymmetry.  For  each 
of  the  resonance  regimes;  analytical  expressions  are  obtained  that  describe 
the  quasi-steady  vehicle  motion,  and  the  results  are  found  to  be  in  good 
agreement  with  computer  solutions  of  the  equations  of  motion.  The  analytical 
approximations  provide  a  simple  tool  for  predicting  reentry  vehicle  dynamic 
behavior  without  requiring  costly  and  time  consuming  machine  computations. 
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center -of-mass  asymmetry 
aerodynamic  axial  force  coefficient 
aerodynamic  roll  damping  derivative 
aerodynamic  pitch,  damping  derivative 
aerodynamic  normal  force  derivative 
aerodynamic  trim  force  coefficient 
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angular  momentum;  altitude 
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pitch  or  yaw  moment  of  inertia 
roll  moment  of  inertia 

static  margin  (distance  of  center  of  pressure  aft  of 
center  of  mass) 
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moment 
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dynamic  pressure 
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time 
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ballistic  coefficient, 
path  angle 
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roll  damping  coefficient, 
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I  U 
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roll  axis 
yaw  axis 

total  angle  of  attack 


roll  moment  coefficient. 


pitch  moment  coefficient. 


<lSc 
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pitch  damping  coefficient, 


s§d 
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pitch  axis 

atmospheric  density 


trim  roll  moment  coefficient, 
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trim  pitch  moment  coefficient, 


roll  orientation  relative  to  wind 
precession  angle 
angular  velocity 

roll 

yaw 

pitch 

oscillatory  component 

time  derivative 
quasi-steady  component 


I.  INTRODUCTION 


It  has  been  known  lor  some  time  that  missiles  or  projectiles  with  sliglit 
configurational  asymmetries  can  exhibit  erratic  roll  behavior  due  to  gyroscopic 
coupling  of  roll  and  pitch. motions  (roll  resonance)  at-or  near  the  natural  pitch 
frequency  of  the  missile  {(Refs.  .1  and  2).  More  recently,  with  the  advent  ,of 
smaller  and  higher  ballistic  coefficient  reentry  vehicles,  similar  erratic 
behavior  has  been  observed  that  can  have  adverse  and  even  catastrophic 
'  effects -on  the  vehicle  motion  during  reentry  (Refs,.  3  and  4). 

Because  of -the  complexity  of  the  equations  of  motion  when  gyroscopic 
coupling  is  involved,  there  :are  few  analytical  treatments  of  roll  resonance, 

.and  much  of  what  is  known  has  been  obtained  from  computer  solutions  of  the 
general  equations  -of  motion.  Such  studies  provide  detailed  Information  on 
the  motions  .and  trajectories  of  .specific  vehicles  but,  unfortunately,  do  not 
lend  themselves  to -identifying  in  a  .scientific  manner  the  pertinent  parameters 
that  -contribute  to  the  roll  resonance  phenomena.  The  purpose  -of -the  present 
study  is  to  identify  analytically  -the  most  important  of these  parameters 
through  ra  {solution  .of  .a  ^simple  set  .of  equations  -embodying  the  .essential 
features  -of  reentry  vehicle  -dynamics.  The  results  so  obtained  provide  a 
simple  tool  for  predicting  the  occurrence  and  .severity  of  reentry  vehicle  roll 
resonance  as  a  function  of  vehicle  -configuration,  .aerodynamic,  and  trajectory 
parameters. 

The  study  treats,  separately,  -three  regimes  of  resonance  Identified  in 
Ref.  4,  in  the  -order  in  which  they  might  be  -encountered  by  ;a  typical  reentry 
vehicle.  First,  the  high  altitude  resonance  when  the  vehicle  roll  .rate  first 
interacts  with  the  critical  roll  rate  is  -examined,  and  the  conditions  necessary 
for  roll  lockin  and  breakout  due  to  a  single  mass  asymmetry  are  described. 
Next,  the  steady  subre  sonant  condition  exhibited  by  a  vehicle  with  a  -single 
mass  asymmetry  when  the  roll  rate  lies  below  the  critical  roll  rate  is 
described  analytically,  and  the  results  are  compared  with -computer  solutions 
of  the  three  -degree -of-freedom  equations  of  motion.  Finally,  an  approximate 


-i- 


method  is  presented  for  predicting  the  transient  roll  spinup  to  resonance 
and  concomitant  angle  -of-attack  divergence  induced,  by  a  compound  asym¬ 
metry,  and  these  results  are  also  compared  with  computer  solutions  of  the 
equations  of  motion. 


4} 


H.  COORDINATES  AND  EQUATIONS  OF  MOTION 


It  is  convenient  for  the  purpose  of  analysis  to  consider  only  three - 
degree -of-freedom  rotational  motion,  since  the  roll  resonance  phenomenon 
is  basically  gyroscopic  in  nature.  There  is,  however,  a  normal  force  con¬ 
tribution  to  pitch  or  yaw  damping  from  lateral  motion  of  the  center  of  mass, 
which  can  be  accounted  for  in  specifying  the  aerodynamic  moments,  so  that 
the  three  moment  equations  effectively  account  for  five -degree -of-freedom 
motion  for  an  axisymmetric  vehicle. 

The  rotational  motion  is  described  in  terms  of  the  Euler  angles  i|j,  <p, 

0,  as  shown  in  Figure  iy  which *4e scribe  the  position  of  the  body-fixed  axes 
x,  y,  z  relative  to  a  set  of  inertial  axes  X,  Y,  Z.  It  is  assumed  that  the  inertial 
frame  of  reference  is  moving  with  the  vehicle  center  of  mass  along  the  velocity 
vector,  which  coincides  with  the  direction  of  the  X-axis,  and  that  the  vehicle 
axis  is  the  body-fixed  axis  x  or  4.  The  angles  i|i ,  <f>,  0,  which  are  rotations 
about  the  axes  X,  x,  £  ,  respectively,  then  represent  precession  about  the 
velocity  vectpr,  roll  orientation  relative  to  the  wind,  and  total  angle  of  attack, 
respectively.  The  mutually  perpendicular'  axes  £,  r|,  £,  which  pre cess  about 
the  velocity  vector,  are  the  pitch,  yaw,  and  roll  axes,  respectively,  where 
the  pitch  plane  is  defined  as  the  plane  of 'total  angle  of  attack.  These  axes 
differ  from  the  usual  body-fixed  axes  x,  y,  %  in  that  the  pitch  axis  £  (the  line 
of  nodes  in  the  Eulerian  system)  is  always  perpendicular  to  the  direction  X 
of  the  velocity  vector,  thereby  defining  the  plane  of  the  total  angle  of  attack, 
whereas  the  body-fixed  axes  y  and  z  rotate  in  this  system  about  the  roll 
axis  (xor  £)  with  angular  rate  <f>. 

If  the  principal  moments  of  inertia  about  the  £,  r\,  £  axes  are  I,  I,  Ix, 
respectively  (pitch,  yaw,  roll),  the  angular  velocities  of  the  £,  r|,  £  axes 
and  the  angular  momentum  about  them  are  (Ref.  5) 


* 


Wg  =  0 

wk  =  ^  sin  0 
*1 

«g  =  4*  cos  0 


hg  =  16 

h  =  ii  sin  0 
T1 


> 


hg  =  1^(0  +  4*  cos  0) 


(1) 


where  <<>g  and  are  the  pitch  and  yaw  rates,  respectively,  and  the  roll  rate 
is  the  sum 


p  =  f  +  +  cos  0 


The  moment  equation  of  motion  is 


M  =  h  +  w  X  h 


which,,  from  Eq.  (1)  with  the  definition,  Eq.  (2),  has  the  components 


1.2 


Mg  =  10  +  Ixp+  sin  0  -  1+  sin  0  cos  0 


(+  sin  0)  +  16+  cos  0  -  Ixp0 


"11  "  dt 

=  Ax  dT 


My  =  I  ^ 


The  moments  Mg,  M^,  Mg  are  the  aerodynamic  pitch,  yaw,  and  roll 
moments,  respectively,  acting  on  the  vehicle. 


(2) 


(3) 


(4) 
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HI.  SINGLE  MASS  ASYMMETRY 


•> 


A.  MOMENT  EQUATIONS 

It  will  be  assumed  initially  that  the  only  asymmetry  that  gives  rise  to 
pitch  roll  coupling,  is.  a. displacement  of  the  vehicle  center  of  mass  from  the 
aids  of  symmetry  (center  of  pressure)  of-the  vehicle  and  that  this  displace¬ 
ment  is  small  enough  that  the  re  suiting  cros  s  products  of  inertia  can  be 
ignored.  Denoting  this  displacement  c,  as.  shown  in  Figure  2,  we  can  write 
the  aerodynamic  moments  in  terms  of  normal  and  axial  force  components  for 
an  otherwise  axially  symmetric  vehicle* 


M^  -  C^-  qSi0  +  C^qSc  cos  </>  + 

,2 


CN  * 

'm  ~ 
q  md 


M  =  C  AqSc  sin  <f>  + 
q  A^  r  u 


CN  1 

r»  Cl 

CV^ 


4»  sin  6 


> 


C(  qSa  p 

M^  =.  -  CN  qS0c  sin  ^  - 


(5) 


J 


1  2 

The  term  Cj^  I/md  in  the  first  two  of  Eqs.  (5)  is  the  usual  normal  force 
damping  termer om  lateral  motion  of  the  vehicle  center  of  mass.  It  can  be 
shown  for  an  axisymmetric  vehicle  that  the  lateral  motions  are  proportional 
to  pitch  and  yaw  rotations  about  a  point  a  distance  I/mi  ahead  of  the  center 
of  mass.  Thus  the  three  resulting  moment  equations  account  for  lateral 
translation  of  the  vehicle  and  effectively  describe  five-degree-of-freedom 
motion. 
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where  C  and  are  pitch  and  roll  damping  derivatives,  respectively, 
and  the  angle  f  describing  the  roll  orientation  relative  to  the  wind  is  zero 
in  the  statically  stable  position  when  the  center  of  pressure  is  on  the  leeward 
side  of  the  center  of  mass,  Negative  <f>  produces  a  positive  roll  torque  and 
positive  ^  produces  a  negative  roll  torque.  It  is  shown  in  Section  IU-B  that 
roll  resonance  is  a  condition  during  which  ^  oscillates  about  a  nonzero  value 
of  <p  producing  a  net  roll  torque.  The  magnitude  and  direction  of  this  torque 
are  strongly  influencedeby  the  magnitude  and  time  rate  of  change  of  the  dynamic 
pressure.  The  oscillation  in  0  about  a  steady  value  or  one  slowly  varying  in 
time  is  defined  here  as  lunar  motion  in  that  a  certain  meridian,  of  the  vehicle 
is  preferentially  facing  into  the  wind.  From  this  definition,  as  is  shown  in 
Section  III-B,  lunar  motion  is  a  necessary  condition  for  resonance. 

Substituting  Eqs.  (5)  in  Eqs.  (4)  and  dividing  the  first  two  equations 
by  I  and  the  third  by  Ix,  we  can  write  the  equations  of  motion 

• 

-\0  +  6  cos  0--  v0 
6  sin  $  -  vty  sin  0 

-«0  sin  0  -  «p  =  -jji- 
where  the  coefficients  are 


=  9  +  |Ap4*  sin  0  -  41  sin  0  cos  0 
=  ('!>  sin  0)  +  0+  cos  0  -  pp0 


The  sign  convention  is  such  that  for  a  statically  stable  vehicle  with  positive 
damping  (C  ,  C  <  0),  all  of  the  coefficients  will  be  positive. 

q  p 

B.  QUASI-STEADY  SOLUTION  FOR  HIGH  ALTITUDE 
ROLL  LOCKIN  AND  BREAKOUT 

It  will  be  assumed  that  the  vehicle  motion  during  reentry,  described 
by  [«Mt),  *(t),  0{t)],  is  of  the  form 

tft),e  (t)]  =  [+(t),^t),?(t)]  +  [4>+(t),*+(t),  e+(t)]  ,  (8) 

where  [»Ji(t),  0(t),  0(t)]  represents  a  quasi-steady  component  that  varies 
relatively  slowly  with  time  {of  the  order  of  the  dynamic  pressure)  and 
f  £j.(t),  0+(t)]  represents  an  oscillation  of. higher  frequency  about  the 
quasi-steady  values.  It  is  of  interest  to  examine  solutions  to  the  equations 
of  motion  for  the  quasi-steady  components  only. 

Making  the  small  angle  approximations  sin  8  =  0,  cos  0  =  i,  neglecting 
pitch  damping  v,  and  ignoring  the  term  6  cos  <p,  which  is  quite  small  except 
at  angles  of.  attack  near  static  trim,  we  can  write  the  first  of  Eqs.  (6)  in  the 
form  _  . 


0  +  (X  +  pp+  -  <i>2)0  =  0  (9) 

This  equation  describes  a  nonlinear  oscillation  in  0  about  a  nonzero  quasi¬ 
steady  value  15",  since  8,  by  definition,  is  always  positive.  The  quasi-steady 
values  iJT  and  p  corresponding  to  0*  are  those  satisfying  the  equation 

X  +  ppI-J2  =  0  ,  (10) 

which  makes  0  in  Eq.  (9)  zero.  Since,  in  general,  p  «  1 ,  the  precession  rate 
+  from  Eq.  (10)  is  approximately  equal  to  sbX.*  ^2,  the  pitch  natural  frequency 
of  the  vehicle.  From  the  definition,  Eq.  (7),  \  '  is  proportional  to  the 
square  root  of  the  dynamic  pressure  and -is  a  relatively  slowly  varying  function 
of  time. 
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Hence,  the  criterion,  Eq.  (15),  with  <p  defined  by  Eq.  (14),  is  a  necessary 
condition  for  resonance  for  a  vehicle  having  a  single  mass  asymmetry.  Re¬ 
writing  «  in  Eq.  (14)  in  terms  of  p  using  the  definitions,  Eqs.  (13)  and  (7), 

c 

we  can  write  sin  <j>  in  the  form 


sin  <j>  = 


(</ C)  Xx/a  -  y 


7% 


( 


±£PCT 


Since  roll  damping  e  is  generally  quite  small,  sin  ^  from  Eq.  (16)  is  approxi¬ 
mately  proportional  to  the  slope  of  the  critical  roll  rate  and  is  inversely 
proportional  to  the  square  of  the  critical  roll  rate  and  angle  of  attack. 

If  the  vehicle  has  an  initial  exoatmospheric  roll  rate,  say  pQ,  then  at 
some  point  in  the  trajectory,  say  tj,  the  vehicle  roll  rate  must  intersect  the 
critical  roll  rate  curve,  as  shown  in  Figure  3. 


Figure  3.  High  Altitude  Roll  Resonance 

t 
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At  this  point  the  vehicle  will  lock  into  resonance,  provided  that  the  criterion, 

Eq.  (15),  is  satisfied  and  that  the  vehicle  has  the  proper  roll  orientation  with 
respect  to  the  wind;  i.  e. ,  the  angle  ^rnust  be  such  that  the  induced  roll  torque 
has  the  same  sense  as  the  critical  roll  acceleration  (e.  g. ,  if  pc  is  positive, 
then  ^  as  defined  in  Figure  2  must  be  negative  to  produce  a  positive  roll  torque). 
Hence,  the  sin  ^  criterion,  Eq.  (IS),  is  a  necessary,  but  not  sufficient,  con¬ 
dition  for  resonance  in  that  lockin  depends  also  on  the  phasing  of  <f>  as  the  roll 
rate  approaches  critical.  •  However,  if  "lockin  does  occurs  then  the  vehicle 
must  break  out  of  resonance  when  sin  ^  is  no  longer  real,  which  will  generally 
occur  at  some  later  point,  such  as  t ^  la  Figure  3,  when  p:  is  large  and  8  has 
damped  to  a  sufficiently  small  value;  The  breakout  point  has  particular 
physical  significance.  Referring  to  Figure  2  and  the  third  of  Eqs.  6,  the  roll 
torque  which  produces  the  roll  acceleration  p  is  proportional  to  the  product 
0  sin  0.  When  0  becomes  sufficiently  small,  the  maximum  possible  roll 
torque,  corresponding  to  $  =  90  deg,  is  insufficient  to  increase  p  at  an 
average  rate  equal  to  the  time  rate  of  change  of  the  critical  roll  rate,  and 
resonance  is  no  longer  possible; 

In  order  to  evaluate  sin  from.Eq.  (16),  we  require  a  knowledge  of  the 
angle -of- attack  convergence  and  the  dynamic  pressure  history,  which 
determines  the  magnitude  and  slope  olthe  critical  roll  rate.  The  dynamic 
pressure  is  approximated  assuming  an  exponential  atmosphere  with  constant 
path  angle,  the  details  of  which;  are  discussed  in  Appendix  A.  The  angle -of- 
attack  convergence  is  obtained  from  a  quasi-steady  solution  of  the  equations 
of  motion,  which  is  summarized  in  Appendix  B.  Using  these  results,  we  can 
now  express  the  sin  ^  criterion,  Eq.  (16),  in.terms  of  the  vehicle  and  tra¬ 
jectory  parameters. 

At  t  =  tj  in  Figure  3,  where  the  .roll  rate  first  intersects  the  critical 

curve,  the  roll  rate  is  approximately  equal  to  the  exoatmospheric  value  p.Q 

and  the  angle  of  attack  has  the  approximate  value  given  by  Eq.  (B-13), 

Appendix  B.  Ignoring  roll  dampiing  e  and  substituting  sin  9^  *  8^  from 

Eq.  (B-13)  with  the  slope. dp  /dt  from  Eq.  (A-16),  Appendix  A,  in  Eq.  (16), 

c 

we  obtain  for  sin  <f>, 
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which,  depending  on  whether  or  not  ^  is  real,  i.  e.  ,  !|sin  -  1 »  is  ^ 
approximate  criterion  for  determining  Toll  locldn.  Assuming  that  lockin 
does  occur,  the  breakout  point  can  also  he  estimated  using  Eq.  (16). 

Equation  (B-15),  Appendix B,  describes  the  angle-of-attack  convergence 
while  the  vehicle  is  in  resonance.  Substituting  this  expression  for  0  in 
Eq.  (16),  using  the  values  for  pc  and  dpc/dtir  om  Appendix  A;  as  before, 
we  obtain  for  the  sin  ^  expression  during  locldn 


sin  tp.  I  - 


1/2 
u  ' 
o 


P*  “  P\  exp  [  (b./2)(p  ~  pp  -  (a,/4)p] 


[B  -  i(p)0 


where 
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and  D(x)  is  Dawsori's  integral  defined  by 


D(x)  &  e”x  f  e"  dt  , 

v 

for  which  tabular  values1  are  available  (see,  for  example.  Ref.  6).  Breakout 
must  occur  when  |sin  p|  from  Eq.  (18)  approaches  unity,  as  shown  in 
Figure  3. 

Figure  4  shows  a  plot  of  sin  p  fromEq;  (18)  for  a  typical  reentry  vehicle, 
for -different  values  of  damping,  measured  from  the  point  at  which  the  roll  rate, 
first  intersects  the  critical  curve,  p.Q  =  pc»  Figure  5  shows  a  comparison  of 
roll  rates  computed  from  a  numerical  integration  of  Eqs.  (6),  with  the 
extimated  breakout  points  based  on  the  sin  p  criterion  of  Eq.  (15),  for  the 
values  plotted  in.  Figure  4.  Also  shown  in  Figure  5  is  a  comparison  of  the 
computed  angle -of- attack  history  with  the  quasi-steady  approximation  of 
Appendix  B  for  C  =  -4.  It  is  seen  that  for  greater  damping,  breakout 
occurs  earlier.  This  is  a  consequence  of  the  increased  angle -of-attack 
convergence,  which  causes  |sin  p|  to  increase  more  rapidly  to  its  limiting 
value  of  unity.  It  is  also  seen  that  the  greater  the  damping,  the  more  accurately 
the  breakout  point  is  predicted.  This  is  a  consequence  of  the  assumption  of 
quasi-steady  motion.  It  was  shown  that  the  vehicle  breaks  out  of  resonance 
when  the  maximum  roll  moment  corresponding  to  p  =  90  deg  is  insufficient 
to  accelerate  the  roll  at  the  critical  roll  acceleration.  Since  the  average 
moment  is  computed  based  on  the  average  moment  arm  c  sin  p  (Figure  2), 
the  moment  will  be  overestimated,  depending  on  the  amplitude  of  the  p  oscil¬ 
lation,  since  the  average  value  of  sin  p  at  <f>  =  90  deg  is  less  than  unity;  i.  e. , 
sin  <f>  <sin  f.  The  lower  the  damping,  the  greater  the  amplitude  of  oscillation 
will  be,  which  accounts  for  the  discrepancies  of  Figure  5. 

C.  SUBRESONANT  CONDITION  ■ 

Following  breakout,  which  is  accompanied  by  a  perturbation  in  the 
vehicle  orientation  and  angular  rates,  and  in  the  absence  of  other  asymmetries. 
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Figure  5.  Comparison  of  Quasi-Steady  Approximation  and  Computer 
Calculation  of  Breakout  Roll  Rate  and  Angle -of- Attack  Convergence 

during  High  Altitude  Resonance 
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TOTAL  ANGLE -OF- ATTACK,  deg 


the  roll  rate  has  been  observed  to  become  almost  constant,  decreasing 
slightly,  and  a  condition  of  steady  lunar  motion  exists  at  very  small  angle 
of  attack.  This  is  a  stable  subresonant  condition  that,  owing  to  its  steady 
nature,  can  be  described  quite  accurately  by  a  quasi-steady  solution  of  the 
equations  of  motion. 

Assuming  that  the  angle  of  attack  and  angular  rates  remain  constant, 
i.  e. ,  0  »  0  and  0*  $  «  n  *  0,  and  the  small  angle  approximations  are  valid, 
we  can  write  the  first1  two  of  Eqs.  (6)  in  the  form 

5  COS  ^  -  {X.  +  HLP'ii  "  =  0 

6  sin  ^  *  0 

where  the  bars,  as  before,  denote  quasi-steady  values.  Since  ^  from 
Eqs.  (21)  must  also  be  slowly  varying,  which  is  the  lunar  motion  condition 
d/dt  ^  »  0  or  h  as'Jj,  Eqs.  (21)  may  be  written 


where  represents  the  almost  constant  value  of  p  after  breakout.  For  p+ 
sufficiently  removed  from  p  ,  as  in  Figure  5,  tan  f  f  «i,  cos  f  *  1,  and 
we  can  write  Eq.  (23)  in  the  form 

T1-  -  <P„/PC>2  •  <24> 

trim 


-18- 


where  9trjm  =  5  A  ic  the  static  or  nonrolling  trim  angle  of  attack  due  to  the 
mass  asymmetry.  Thus,  it  is  seen  that  this  trim  is  amplified  by  the  factor 
[l  -  (p*/pc)^j  "*»  which  is  characteristic  of  a  near-resonant  condition.  It 
is  shown  in  Section  IV -B  that  an  aerodynamic  trim  in  conjunction  with  the 
mass  asymmetry,  which  can  cause  spinup  into  low  altitude  resonance,  is 
amplified  by  this  same  factor  during  spinup. 

The  slight  roll  deceleration  during  the  subresonant  condition  may  be 
estimated  from  the  third  of  Eqs.  (6)  upon  substitution  of  sin  f  *  tan  ^  from 
Eq.  (22)  and  IT  from  Eq.  (23),  which  gives 


This  deceleration  is  found  to  be  quite  small  except  as  p^  —  p  . 

Figure  6  shows  a  comparison  of  the  quasi-steady  angle  of  attack  from 

Eq.  (24)  with  that  obtained  from  a  computer  solution  of  the  equations  of 

motion,  for  the  case  C.  =  -4  discussed  previously.  Also  shown  in  the 

mq  <■ 

figure  are  the  computed  roll  rate  and  quasi-steady  approximation  bat  .id  on 
the  breakout  roll  rate  from  Figure  5  with  a  numerical  integration  of 
Eq.  (25). 
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Figure  6.  Comparison  of  Quasi-Steady  Approximation  and  Computer 
Calculation  of  Angle -of -Attack  and  Roll  Rate  during  Subresonance 
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ROLL  RATE,  rad/ sec 


.  IV.  COMPOUND  ASYMMETRY 

A.  MOMENT  EQUATIONS  WITH  AERODYNAMIC  TRIM 

It  has  been  shown  (Ref.  4)  that  if  a  relatively  small  aerodynamic  trim 
force  develops  at  lower  altitude  that  acts  out  of  the  plane  of  the  center-of- 
mass  asymmetry,  a  roll  acceleration  will  be  produced  that  can  be  sufficient 
to  spin  the  vehicle  into  low  altitude  resonance.  This  trim  can  be  represented 
by  the  normal  force  Ft  oriented  at  an  angle  with  respect  to  the  plane  of 
the  center -of -mass  asymmetry,  as  shown  in  Figure  7,  and  acting  at  a  dis¬ 
tance  L  aft  (or  fore)  of  the  center  of  mass.  Referring  to  Figure  2  and 

y  j 

defining  a  trim  force  coefficient  C^  =  F^/qS,  we  write  the  resulting  pitch, 
yaw,  and  roll  moments  in  the  form 

M^  =  CtqSLtcos  (0  +  0Q) 

Mq  =  C^SI^sin  (<f>  +  0Q)  (26) 

=  -CtqSc  sin 

When  these  terms  are  added  to  the  moments  of  Eq.  (5),  the  equations  of 
motion,  Eqs.  (6),  become 

t  cos  (£  +  4>Q)  -  \0  +  6  cos  ^  -  t>0  =  0  +,pp<i»  sin  6 
t  sin  (0  +  0Q)  -  sin  0  =  ^ (^  si11  ®)  +  ©4*  cos  0 
-  cr  sin  <f>Q  -  *0  sin  <f>  -  ep.  =  ^ 
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where  t  and  r  are  defined  by 


C.qSL  C  qSc 

Y  V  ■  ,  .  (28) 

x, 

When  fQ  »  ±90  deg,  a  relatively  small  normal  trim  force  acting  at  a  sufficient 
distance  from  the  center  of  mass  (at  the  fore  or  aft  ends  of  the  vehicle, 
for  example)  will  cause  the  roll  to  accelerate  from  its  relatively  constant 
subresonant  value  to  the  critical  value,  as  shown  in  Figure  8.  On  reaching 
the  critical  roll  rate  the  vehicle  wiil  once  again  lock  into  resonance  with  a 
subsequent  divergence  in  angle  of  attack  and  increase  in  lateral  loading.  The 
same  trim  introduced  at  higher  altitude  prior  to  breakout  from  resonance 
would  cause  the  vehicle,  to  remain  locked  into  resonance. 


Figure  8.  Spinup  to  Low  Altitude  Resonance  from 
Compound  Asymmetry 


The  roll  and  angle -of- attack  behavior  during  spinup  and  resonance 
can  be  obtained  from  an  approximate  quasi-steady  solution  of  Eqs.  (27), 
similar  to  that  used  in  treating  the  high-altitude  resonance  problem-.  It  is 
convenient  to  treat  the  spinup  and  resonance  phases  separately. 

B.  SPINUP  TO  RESONANCE 

The  mechanism  for  spinup  to  low  altitude  resonance  with  the  presence 

of  a  body-fixed  trim  is  similar  to  the  roll  acceleration  during  high  altitude 

resonance  due  to  a  single  mass  asymmetry.  In  both  cases,  the  roll  torque 

is  the  normal  force  from  angle  of  attack  acting  on  the  center-of-mass 

asymmetry  as  its  moment  arm.  During  high  altitude  resonance,  the  angle 

of  attack  is  the  residual  from  the  eroatmospheric  reentry  value,  which  has 

not  completely  converged,  whereas  during  spinup,  the  angle  of  attack  is 

induced  by  the  aerodynamic  trim.  In  order  that  the  roll  torque  be  sustained, 

the  vehicle  must  have  a  preferential  roll  orientation  relative  to  the  wind, 

which  is  the  lunar  motion  condition  0  =  0. 

Again,  if  we  consider  only  the  quasi-steady  vehicle  motion  by  neglecting 

•  •• 

the  angle -of-attack  oscillations  exemplified  by  0  and  0,  and  make  the  small 

angle  approximations  as  before,  the  lunar  motion  condition  requires  that 
•  _ 

4»  »  p,  and  we  can  write  Eqs.  (27)  as 

t  cos  (0  +  0q)  -  X.0  +  6  cos  ^  =  -p^(l 

r 

r  sin  -  vp  0  =  0p. 

♦ 

-or  sin  -  «0  sin  ^  =  p 

We  can  solve  these  equations  approximately  for  the  quasi-steady  values  p 
and  0,  for  the  case  =  ±90  deg,  by  neglecting  the  terms  6  and  < r,  which  are, 
in  general,  small  relative  to  other  terms  in  the  respective  equations.  It  is 
shown  in  Section  IV-C,  by  a  comparison  of  the  results  so  obtained  i -1th  a 
numerical  integration  of  Eqs.  (27),  thatthis  is  a  good  assumption  fir  typical 


vehicle  and  trajectory  parameters,,  Using  the  definition  of  p  ,  Eq;  (13), 
we  can  then  reduce  Eqs.  (29)  to 

-r  sin  0  ='\|i  -  (p/pc)2Je 

T  COS  0  =  0  p  j 

♦  ' 

-  «0  sin  0  =  p  „ 

By ^liminating  (F  and  0  between  these  equations  using  the  identity 
sin  0  +  cos  0=  1,  we  obtain  the  nonlinear  differential  equation  for  p 


pV[‘-s/p/]  *ipiF/p/)  '  (31) 

♦ 

This  is  a  cubic  equation  in  p with  the  second  order  term  missing,  which  has 
the  exact  solution 


where 


3  VJkt 


f  -  <p/p7J 


For  typical  vehicle  and  trajectory  parameters,  X  «  1  for  |p/p.  |  <0.  9,  and 
Eq.  (32)  can  be  expanded  in  powers  of  X,  which  yields  the  simple  result 


-  <F/p/] 
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♦ 


This  indicates  that  the  second  term  in  Eq.  (31.)  can  he  neglected  for 
j  p/p  |  <  o.  9.  Since  p^  and  the  coefficients  «,  T»  and  X  are  proportional  to 

the  dynamic  pressure,  which  is  lime  dependent,  Eq.  (34)  can  he  solved  for 
p  by  a-simple  numerical  integration.  Substituting  Eq.  (34)  in  Eqs.  (30.).,, 
we  find  the  corresponding  angle -of-attack  behavior  during  spinup  to  be 

(35) 

This  is  the  same  trim  amplification  factor  derived  in  Section  IH-C  for  the 
subresonance  condition  with  a  single  mass  asymmetry. 

C.  LOW  ALTITUDE  RESONANCE  AND  ANGL.E-OF-  ^ 

ATTACK  DIVERGENCE 

As  the  roll  rate  approaches  the  critical  value  Eqs.  (34)  and  (35)  no 
longer  apply,  and  it  is  essential  to  include  pitch. damping  v  to  obtain  the  roll 
behavior  during  resonance;  As  with  the  usual  resonance  phenomenon,  the 
angle  of  attack  exhibits  a  singularity  in  the  al  sence  of  damping  as  p  — plc* 
Therefore.,  retaining  v,  we  may  write  Eqs.  (29,)'with'0o  -  ‘90  deg  in  the  form 


6  cos  -0—  t  sin  0  ~  X 


t  cos  0  =  (vp  4-  p  )  '6 
—  <  tr  —  «c0  sin  0  =  p 


<36) 


As  p  -~p  ,  the  first  of  Eqs.  (36)  gives  tan  0  —  6/t  =  (cA/Ct)(c//b*t),  which, 
in  general ,  is  a  very  small  quantity  for  trim  moments  (C^L^)  of  sufficient 
magnitude  to  spin  the  vehicle  into  resonance  before  impact.  Consequently, 
sin  0—0  and  cos~0  —1,  which  gives  for  9  during  resonance,  from  the 
second  of  Eqs.  (36), 
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Noting  that  the  nonrolling  trim  angle  of  attack  0^  due  to  the  trim  t  is  t/X, 
we  find  that  the  amplification  of  the  nonrolling  trun  during  resonance  is 


i.  (1  ~  rtfl 

et  =  vpc  +  pc 


(38) 


When  p  =  0,  at  peak  dynamic  pressure,  this  expression  reduces  to 
c 


r<i-ycvi‘/2 

2  41 

L  J 

[(CN<(I/mdi)]dJ 

(39) 


which,  .with  the  exception  of  I  in  place  of  I  -  Ix  in  the  denominator,  is  the 
identical  expression  derived  in  Ref.  4  using  linear  theory.  The  trim 
amplification  factor,  Eq.  (38),  includes  the  dependence  of  trim  amplification 
on  the  slope  of  the  critical  roll  rate. 

The"  foregoing  results  provide  a  simple  means  for  predicting  the  roll 
acceleration  and  angle -of-attack  divergence  as  a  function  of  time  (or  altitude) 
after  the  introduction  of  an  aerodynamic  trim.  Because  of  the  dependence 
of  the  angle  of  attack  during  resonance  on  the  slope  of  the  critical  roll  rate 
curve  and  on  pitch  damping  (which  is  proportionalto  vehicle  velocity),  the 
maximum  angle  of  attack  will  depend  on  the  altitude  at  which  the  vehicle 
spins  into  resonance.  The  approximate  angle -of- attack  behavior  can  be 
constructed  from  the  results  of  Eqs.  (35)  and  (37). 

Figure  9  shows  a  typical  comparison  of  the  quasi-steady  approxima¬ 
tion  with  a  computer  calculation  of  the  roll  acceleration  to  resonance  after 
the  introduction  of  a  body-fixed  aerodynamic  trim  90  deg  out  of  the  plane  of 
the  center-of-mass  asymmetry.  Figure  10  shows  a  similar  comparison  of 
the  corresponding  angle -of-attack  divergence. 
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Figure  9.  Comparisonof  Quasi-Steady  Approximaticn  with  Computer 
Calculation  of  Roll  Rate  Transient  after  Introduction 
of  Aerodynamic  Trim 
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Figure  10;  Comparison  of  Quasi-Steady  Approximation  with  Computer 
Solution  of  Angle -of -Attack  Divergence  after  Introduction 
of  Aerodynamic  Trim 
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V.  CONCLUSIONS 


It  has  been  demonstrated  that  reentry  vehicle  roll  resonance  can  be 
described  by  a  solution  of  the  moment  equations  of  motion  for  the  quasi-steady 
or  average  values  of  the  vehicle  orientation  and  angular  rates,  without  regard 
for  perturbations  in  the  motion  about  these  values.  This  results  in  a  con¬ 
siderable  simplification  of  the  equations  of  motion. 

The  analysis  has  been  applied  to  different  regimes  of  reentry  vehicle 
roll  resonance  behavior,  previously  identified  through  computer  solutions 
of  the  general  equations  of  motion  for  particular  vehicles  and  trajectories. 

The  analytical  results  provide  a  simple  tool  for  studying,  parametrically, 
reentry  vehicle  motion  for  different  vehicle  configurations  and  trajectories, 
and  they  serve  to  identify  the  important  parameters  and  their  relative 
influences  on  roll  resonance  phenomena. 

It  should  be  possible  to  extend  the  analysis  to  include  other  aerodynamic 
forces,  such  as  pure  roll  torques  or  forces  due  to  ablation,  or  mechanically 
induced  forces  that  might  be  desirable  for  control  of  the  vehicle  motion. 


* 
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APPENDIX  A 


\J*  - 


CRITICAL  ROLL  RATE  APPROXIMATION  FOR 
EXPONENTIAL-ATMOSPHERE 


2 

Using  the  approach  of  Allen  and  Eggers,  we  can  relate  the  critical 
roll 'rate  and  its  derivatives  to  vehicle  and  trajectory  parameters  as  follows. 
It  is  assumed  that  atmospheric  density  p  varies  exponentially  with  altitude  h 
according  to 


P  = 


(A-l) 


where  p-  and  H  are  constants.  Neglecting  acceleration  from  gravity  and 
assuming  a  constant  axial  force  coefficient  C^,  we  can  write  the  vehicle 
deceleration  from  drag  in  the  form 


(A -2) 


where  |3  is  the  ballistic  coefficient  defined  by 


=  J2£L 

-cAs 


For  a  constant  path  angle  y,  the  vertical  descent  rate  is 


dh 

dt 


=  u  sin  6 


(A-3) 


(A -4) 


7 - 

J /.  H.  Allen  and  A.  , 3".  Eggers,  Jr. ,  A  Study  of  The  Motion  and  Aerodynamic 
Heating  of  Ballistic  Missiles  Entering  the  Earth's  Atmosphere  at  High  Super¬ 
sonic  Speeds,  NACA  Report  No.  1381  (1958). 


Differentiating  Eq.  (A-l)  with  respect  to  time  and  making  use  of  Eq.  (A-4), 
we  find  that  the  rate  at  which  density  increases  during  descent  is 

d£  pu  gin  i  (A  5) 

dt  H 

Eliminating  t  between  Eqs.  (A -2)  and  (A-5),  we  find  that  the  velocity  change 
with  density,  in  terms  of  the  trajectory  parsuneter 

*  ^irfib  •  ,A-6) 


is 


du  _  ua 
dp  -  "2”  ' 


which  may  be  integrated  to  give 


«(a/2)p 

u  =  uQe  '  •  , 


where  uq  is  the  reehtry  velocity. 

The  dynamic  pressure  as  a  function  of  density  is  then 


1  2  i  2  -ap 

q  =  2  Pu  =2  puoe 


which  has  a  maximum  value 


(A-7) 


(A-8) 


(A-9) 


2 


u 


‘q*  =  27£ 


> 


A-2 


(A-10) 


at  an  altitude 


h^Hlnp^/p,,, 


(A- 


where  the  density  is 


(A- 


and  the  velocity  is 


u*  =  uq/n/« 


(A- 


The  critical  roll  rate,  p  ,  defined  by 

C 


1./2 


(A- 


may  now  be  expressed  as  a  function  of  atmospheric  density  and  trajectory 
parameters  by  substituting  for  q  from  Eq.  (A-9),  which  gives 


Sf  v  1./2 


■K) 


%  i/2e-(a/2)p 

Nf2 


(A- 


The  time  rate  of  change  of  pc»  which  is  also  of  interest,  "is  found  to  be 


The  atmospheric  density  is  related  to  altitude  through  Eq.  (A-l)  and  may 
be  expressed  as  a  function  of  time  through  the  expression 


dt  H 


^  =  ^2.e^a/2)Psinv  , 


which  may  be  integrated  to  give 


t  -  t  = 


H 


o  u  sm  y 
o 


’  p  .  <a./2)Cp  .  pc)  <a/2)2(p2  -  p2) 
ln  p  +  1  *  i !  +  2  •  Tv 


(A-17) 


*  •  « 


(A-18) 


APPENDIX  B 


ANGLE-OF-ATTACK  CONVERGENCE  DURING  ROLL  LOCKIN 


We  can  estimate  the  angle -of -attack  convergence  through  a  quasi¬ 
steady  solution  of  the  equations  of  motion, 

•  •  *  •  O 

-  \0  +  6  cos  Q  -  v0  =  (§  +  pp4*  sin  0-4*  sin  0  cos  0 
6  sin  ^  -  vty  sin  0  =  t^(4  sin  ©)■+  ©4*  cos  0  -  jxp0 

-  *0  sin  ^  -  cp  = 

using  the  exponential  atmospheric  density  approximation  discussed  in 
Appendix  A.  During  resonance,  the  vehicle  is  in  a  state  of  lunar  motion, and, 
for  small  0,  the  quasi-steady  or  average  values  of  p  and  4*  are  equal  to  the 
critical  roll  rate;  i.  e. , 


(B-i) 

» 


(B-2) 


For  a  vehicle  with  a  large  pitch-to-roll  inertia  ratio,  p  is  small  compared 
with  cos  0,  and  the  term  pp0  in  the  second  of  Eqs.  (A-i)  can  be  neglected 
compared  with  04J  cos  0.  We  can  then  write  this  equation,  with  the  condition 
of  Eq.  (B-2),  in  the  form 


(B-3) 


B-i 
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If  we  substitute  the  expression  for  sin  ^  from  the  third  of  Eqs.  (B-l)  into 
Eq.  (B-3),  with  the  small  6  approximation  and  the  conditions  Eq.  (B-2), 
Eq.  (B-3)  becomes 
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.  2 
sxn 


*) 


2— 

+  vp.  sin  9  = 
*c 


(B-4) 
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This  equation  can  now  be  integrated  for  pcsin  9  and  has  the  familiaT 
solution 


2—  , 
pcsin  9  =  exp  | 


[\^  v(t'>  dt]  (Pcate20)tn1  -  4  f  «?[/  W‘”)  dt] 


where  time  t  =  tj  corresponds  to  the  point  at  which  the  vehicle  first  loclcs 

into  resonance,  i.  e. ,  when  the  roll  rate  p  intersects  the  critical  curve  p  . 

c 

The  angle  of  attack  9  at  this  point  can  be’  estimated  from  the  exoatmospheric 
vehicle  motion  as  follows. 

Before  reentry,  in  the  absence  of  external  moments  and  pitch  accelera¬ 
tion  9,  the  vehicle  will  have  some  exoatmospheric  roll  and  precession, motions 
pQ  and  4*o  at  angle  of  attack  9q,  which  are  gyroscopically  coupled  according  to 


\ 


wo 
cos  9 


o 


(B-6) 


Since  4*  and  pQ  will,  in  general,  be  unequal,  the  position  of  the  windward 
meridian  £  will  rotate  continuously  about  the  vehicle,  and  a  steady  roll  torque 
cauiot  persist  as  it  would  in  the  case  of  lunar  motion.  However,  there  can  be 
a  net  roll  torque  depending  on  the  phasing  of  (f>  relative  to  angle -of-attack 
oscillations  9.  Nevertheless,  in  view  of  the  small  magnitude  of  the  aerodynamic 
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forces  at  these  altitudes,  It  is  assumed  that  p.  =  pQ  -  constant  until  the  roll 
’  rate  intersects  the  critical  curve.  If  we  neglect  roll  damping  the  third  of 

Eqs.  (B-l)  requires  that  sin  f  =  constant.  Since,  from  considerations  of 
*•  the  quasi-steady  solution  of  the  pitch  moment  equation  discussed  in  the  body 

of  this  report,  the  quasi-steady  precession  rate  4*  will  be  equal  to  the  natural 

1/2  — 
pitch  frequency  X  ,  which  is  very  nearly  p  ,  Eq.  (B-3)  with  sin  ^  =  0 

describes  the  angle -of- attach  convergence  prior  to  high  altitude  resonance. 

Integrating  this  equation  from  0  to  t,  we  obtain 

0  <  t<  t^ 

(B-7) 


9 _  ;  2 

p„  sin  0  =  ^  sin  0  exp 
rc  o  o 


[-/twt',dt'] 


We  can  evaluate  the  integrals  of  Eqs.  (B-5)  and  (B-7)  using  the  expo¬ 
nential  atmosphere  approximations  for  q  andu  derived  in  Appendix  A.  Changing 
the  independent  variable  from  t  to  p,  using  the  relation 


-(a/2)p 


smy 


(B-8) 


derived  in  Appendix  A,  and  neglecting  roll  damping  in  Eq.  (B-5),  we  find 
that  the  integrals  are 


(B-9) 
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where 
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and  D(x)  is  Dawson! s  integral  defined  by 


2  /.x  2 

D(x)  =  e  x  /  e 


dt 


(B-12) 
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for  which  tabular 'values  are  available.  The  angle-  of-attack  expressions, 
Eqs.  (B-5)  and  (B-7),  can  now  be  written  in  terms  of  these  results.  Of 
particular  interest  is  the  angle  of  attack  when  the  roll  rate  first  intersects 
the  critical  curve  at  t  =  t^.  At  this  point  the  critical  roll  rate  will  equal 
the  exoatmospheric  value  pQ,  which  was  assumed  to  remain  constant  in  the 
interval  Os  t  s  t^.  Therefore,  from  Eq.  (B-7)  at  t  =  t^,  using  the  relation, 
Eq.  :(6),  wei'obtain 


_Ji _ 

cos  9 

o 


sin  6q*  exp  [-  (b/2)(pi 


(B-13) 
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See,  for  example,  Handbook  of  Mathematical  Functions  with  Formulas, 
Graphs  and  Mathematical  Tables',  Applied ^Mathematics  Series  55,  National 
Bureau  of  Standards  (1964). 
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